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ABSTRACT 

We perform a statistical analysis of the peculiar velocity field around dark mat- 
ter haloes in numerical simulations. We examine different properties of the infall of 
material onto haloes and its relation to central halo shapes and the shape of the large 
scale surrounding regions (LSSR). 

We find that the amplitude of the infall velocity field along the halo shape minor 
axis is larger than that along the major axis. This is consistent for general triaxial 
haloes, and for both prolate and oblate systems. We also report a strong anisotropy 
of the velocity field along the principal axes of the LSSR. 

The infall velocity held around dark matter haloes reaches a maximum value along 
the direction of the minor axis of the LSSR, whereas along the direction of its major 
axis, it exhibits the smallest velocities. We also analyse the dependence of the matter 
velocity field on the local environment. The amplitude of the infall velocity at high 
local density regions is larger than at low local density regions. The velocity field tends 
to be more laminar along the direction towards the minor axis of the LSSR, where the 
mean ratio between flow velocity and velocity dispersion is of order unity and nearly 
constant up to scales of 15 Mpc h _1 . 

We also detect anisotropies in the outflowing component of the velocity field, 
showing a maximum amplitude along the surrounding LSSR major axis. 
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1 INTRODUCTION 

Given the current paradigm for structure formation in the 
Universe, mass is accreted onto haloes from the network 
of filament walls and voids. Therefore, a strong correla- 
tion between halo properties and large scale dynamics is 
expected in this scenario. Even though statistical isotropy 
at very large scales is a fundamental assumption in cosmol- 
ogy widely verified in galaxy surveys, on smaller scales of a 
few Mpc, the structure of the mass distribution departs sig- 
nificantly from spherical symmetry with the usual presence 
of walls, filaments and voids. This anisotropy in the mat- 
ter distribution has been extensively studied, as it is related 
to the search for large scale structures in the environment 
close to galaxies. For exam ple, both observational studies 
and numerical analysis (e.g . West 1994; Plionis fc Basilakosl 



I2OO2I ; iKitzbichler fc Saurerj 120031 ; iFaltenbacher et alj l2002h 

reported that galaxies tend to be aligned with their neigh- 
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bours supporting the vision of anisotropic mergers along fil- 
amentary structures. When considering preferential direc- 
tions along the large scale structures, it is natural to un- 
derpin long filamentary structures connecting large clusters. 
Regarding the evolution of structure in simulations, the flow 
of particles and haloes within these filaments is apparent and 
can originate preferential directions in the ve locity field (e.g. 
iGonzalez fc Padillall20icl : iPivato et al.ll2006l ). 

On the other hand, the anisotropic distribution of struc- 
ture around haloes, induces the emergence of vortical mo- 
tions on the halo velocity field, through tidal interactions. 
Previous works have shown that the distributions of spin 
vectors are not random. For example, the spin of haloes in 
simulations t end to point in t he direction orthogonal to the 
filaments (see lPaz et alj2008l and references therein). Com- 
bined with the results su ggesting that halo spins are more 
sensiti ve to recent infall (|van Haarlem fc van de Weygaert] 
Il993al ). these alignment properties fit well with accretion 
along preferred directions. In addition, several authors have 
shown that halo shapes are aligned with the surrounding 
large scale matter distribution, with increasing ellipticity as 
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the mass increases (see for instance IPaz et aTll200d ). There- 
fore halo shapes are likely to depend on the accretion process 
driven by halo mergers. Following this line, it is expected a 
strong correlation between the velocity field in the neighbor- 
hood of haloes and the directions defined by the halo shape 
axes. 

Most of these previous studies focused on the fact that 
alignments and preferential directions are consequences of 
the formation process of haloes. However, the effects of such 
preferential directions on the peculiar velocity field at large 
scales have not been so extensively addressed. In a hierar- 
chical clustering scenario these peculiar velocities are gener- 
ated by inhomogeneities in the distribution of matter in the 
universe. The nature of this velocity field depends on the 
local density. High density regions show random motions 
typical of virialised objects whereas low density environ- 
ments, on the other hand, are more likely to show streaming 
motions: objects falling towards larger potential wells con- 
stantly increasing the am plitude of their clustering strength 
l|Diaferio fc Gerlerlll997f ). 

In this work we use a numerical simulation of the stan- 
dard scenario of Cold Dark-Matter with a cosmological con- 
stant (ACDM) to explore different features of the peculiar 
velocity field and its relation to haloes and their large-scale 
surrounding region (referred to as LSSR hereafter). In Sec- 
tion 2 we describe the analysis of the simulation, Section 3 
deals with the anisotropics of the velocity field in relation 
to halo and LSSR shape and orientation. We also consider 
in this section the infall of particles and haloes in different 
local density regions. In Section 4 we explore the turbulence 
of the flow along different directions and finally in Section 
5 we analyse the outflow of particles from haloes and its re- 
lation to the LSSR. We present our conclusions in Section 
6. 



2 DATA 

Throughout this work we use a collisionless numerical sim- 
ulation of 512 3 particles covering a periodic volume of 500 3 
(h _1 Mpc) 3 . We assume a spatially flat low-density Universe, 
with a matter density fi m = I — £7a = 0.28, Hubble con- 
stant H = 74 km s _1 Mpc -1 , and normalisation parameter 
(78 = 0.8. With these parameters the resulting particle res- 
olution is m p — 7.2 x IO lo M0h _1 . The initial conditions 
at r edshift ~ 50 were generated with the GRAFIC2 pack- 
age (|Bertsc hingcr 2001), which also computes th e transfer 
function as described in lMa fc Bertschinger] (jl995j) . The sim- 
ulation run was performed using the secon d version of the 
GADGET code developed bv lSpringel (|2005T ). wi th a gravita- 
tional softening of 0.03h _1 Mpc chosen following lPower et al.l 
|2003l ). 

The identification of particle clumps to be associated to 
dark matter haloes was carried out by means of a standard 
friends-of-friends algorithm with a percolation length given 
by I = 0.17 u~ 1/3 , where v is the mean number density of 
dark matter particles. We identified approximately 470, 000 
dark-matter haloes with at least 10 particle members, re- 
sulting in a minimum halo mass of M m i n = 7.2 IO^MqIi -1 . 
It should be noticed that this algorithm does not resolve 
substructures within the dark-matter haloes. 

From this point on, we will centre most of our analy- 



sis on samples of central haloes with different masses. We 
draw five halo subsamples with equal number of haloes of 
different mass starting at a minimum of 50 particles or 
3.5 x lO 12 M h _1 (the lowest mass sample considered in- 
cludes haloes of up to 64 particles each). The subsample 
with the most massive haloes starts at 260 particle mem- 
bers corresponding to 1.9 x lO 13 M0h _1 . 

Throughout this work, we estimate statistical errors by 
means of the jacknife technique, which has been shown to 
provide equivalent results to error estimates obtained from 
the variance of a large number of i ndependent simulation s 
(see for instance. [Croton et alj2004lPadilla fc Baughll2003h . 
These results indicate that jacknife errors provide a reason- 
able estimate of statistical uncertainties and cosmic vari- 



3 VELOCITY FIELD ANISOTROPIES 

In this section, we analyse anisotropies of the radial velocity 
field relative to the centre halo shape axes. We also anal- 
yse the dependence of this velocity field with respect to the 
surrounding mass distribution at small and large scales. 



3.1 Dependence on the halo shape 

Several authors JWarren et all [l992l; iThomas et al.1 Il998l; 



Hopkins et al 
20051 : 



2005; 



Allgood et alJ l200d; iKasun fc Evrardl 
Lau et al.l 120101 : IPaz et al.l 120061 . and references 



therein) have analysed the properties of halo shapes using 
the best-fitting ellipsoid to the spatial distribution of the 
halo material. Following this standard method, we calcu- 
late the shape tensor for each dark-matter halo using the 
positions of its particle members. This can be written as a 
symmetric matrix, 



(l/Nh 



E 



(i) 



where X a i is the i* component of the displacement vec- 
tor of a particle a relative to the centre of mass, and Nh 
is the number of particles in the halo. The matrix eigen- 
values correspond to the square of the axis (a, 6, c were 
a > b > c) of the characteristic ellipsoid that best de- 
scribes the spatial distribution of the halo members. The 
principal axes of this ellipsoid are represented by the eigen- 
vectors a, b and (Q. Haloes with b/a < c/b are associated 
to general triaxial ellipsoids with prolate shape, while the 
opposite case, b/a > c/b, corresponds to predominantly 
oblate systems. Numerical results based on this shape analy- 
sis have shown that haloes exhibit aspherical shapes, with a 
slight preference toward prolate forms. Moreover, it has been 
shown that the ir orientations are related to their surround- 
ing structures ((Colberg e t al. 2005; Kasun & Evrard 20051 ; 
| Basilakos et aljhood: lAllgood et al.ll200d; [Xitav et al.ll200d; 
Aragon-C alvo et al.1 120071 ; iBrunino et al.l 120071 ; iBett et al.l 
20071 ; IZhang et al.1 120091 '). Several authors have found that 



1 A fixed b/a = 1 with an arbitrary value of c/b, corresponds to 
perfect oblate ellipsoids. On the other hand, systems with fixed 
c/b = 1 are perfect prolate ellipsoids 
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Figure 1. Upper panel: histograms showing the shape (left) 
and eccentricity distribution (right) of haloes (red) and the LSSR 
(black). Lower panel: normalised distribution of the angle be- 
tween haloes and LSSR minor axes (black). The corresponding 
distribution for major axes is shown in red lines. 



dark matter haloes tend to be more prola te and aspherical 
when larger halo mas ses are considered l|Kasun fc Evrardl 
120051 : IPaz et alJ 120061 . and references therein). This trend 
could be easily understood as the result of an accretion pro- 
cess driven by halo mergers. Therefore, we have a physical 
motivation to expect a connection between the halo sur- 
rounding velocity field and the particular directions defined 
by the halo shape axes. 

The resulting halo ellipticities, which are defined as the 
ratio c/a, are shown in red line s in th e right upper panel 
of Figured] Following IPaz et alJ (2006) we separate prolate 
and oblate ellipsoids through a triaxiality parameter "P" 
defined as: 



P = log{{c/b)/{b/a)). 



(2) 



where systems characterised by P > (P < 0) correspond to 
halo with prolate (oblate) shapes. In the left upper panel of 
Figure [1] we show the distribution of triaxiality parameter 
values (red line) derived for haloes. As can be noticed in 
this figure haloes show a trend towards prolate shapes, in 
agree ment with previous results (see for instance I Paz et all 

l2ooeh 

In order to characterise velocity field anisotropies 
around haloes, we calculate the mean radial velocity of par- 
ticles taking into account their positions relative to the dark 
matter halo major axes. We compute radial velocities rel- 
ative to the halo centre, with positive velocities indicat- 
ing mean infall flux and negative values, outflow velocities. 
Along this subsection we analyse elongated haloes, charac- 
terised by small ellipticity parameters (e = c/a < 0.5). From 
this subsample we select either prolate or oblate haloes, cor- 
responding to triaxiality parameters P > 0.1 or P < —0.1. 
From this point on we will simply refer to these two halo 
subsamples as prolate and oblate haloes. 



For prolate haloes, we define 9 a as the angle subtended 
by the halo-particle relative position and the major axis. We 
calculate the average radial velocity as a function of halo- 
particle separation in three cases; i) for particles close to 
the major axis (8 a < 7r/6), ii) for particles roughly perpen- 
dicular to the major axis (tt/3 < 9 a < w/2) and hi) for 
particles at intermediate positions (tt/6 < 8 a < t/3). The 
resulting velocity curves are shown in the left panels in Fig- 
ure [21 where the upper panel corresponds to low mass haloes 
(rh halo = 4xlO 12 M h- 1 ; 7.2 x lO^M©!!- 1 < m ha i < 4.5 x 
lO 12 M0h _1 ) and the lower panel corresponds to high mass 
haloes (m halo = 6x lO^Moh" 1 ; 1.9x IO^MqIT 1 < m halo ). 
Red lines in the left panels Figure [2] show the velocity field 
traced by particles along the major axis ; cyan lines corre- 
spond to particles in intermediate positions and black lines 
correspond to velocities of particles with positions perpen- 
dicular to the major axis. Similar results were obtained for 
the intermediate mass bins. 

In a similar fashion, we derive the velocity field for par- 
ticles along the minor axis and in the perpendicular direc- 
tion for oblate haloes. We define 9 C as the angle subtended 
by the halo-particle relative position and the minor axis 
direction. The mean velocities obtained are shown in the 
right panels of Figure [5] where black lines show the velocity 
field traced by particles along the minor axis (8 C < tt/6), 
red lines correspond to velocities of particles with positions 
perpendicular to the minor axis (tt/3< 8 c < tt/2), and 
cyan lines correspond to particles at intermediate positions 
(n/6< 8 C < tt/3). We use the same halo mass ranges in both 
left and right panels of Figure 

As can be seen in Figure [2] velocities are positive over 
the entire scale range, indicating that the surrounding mate- 
rial is infalling onto haloes regardless its mass or direction. 
It is worth mentioning that the larger infall velocities are 
obtained for the more massive haloes , in agreement with 
previous results ijCeccarelli et al.ll2005h and t heoretical pre- 
dictions l|Pivato et al.ll2006l : ICroft et al ]| 19991 . and references 
therein) . 

Regarding the analysis of prolate haloes, it can be seen 
in the left panels of Figure [2j that the mass infall ampli- 
tude is sightly smaller in the direction along the major axis 
(red lines) than the corresponding infall obtained over the 
perpendicular plane to major axis (black lines). For oblate 
haloes (right panel on Figure [5J, the largest infall velocities 
are obtained at directions parallel to the minor halo axis 
(black line), which exceed by approximately 40 km/s the 
corresponding velocities obtained along the perpendicular 
direction (red line). Both observed anisotropies for prolate 
and oblate halo subsamples, can be understood as a ten- 
dency of the surrounding mass to fall faster from minor axis 
directions, which are the directions set by the more flattened 
sides of the halo shape ellipsoid. This difference is detected 
up to distances to the halo centre of the order of a few times 
the virial radius (approximately up to distances around 2 to 
3 Mpc h" 1 ). 

3.2 Dependence with the surrounding structure 

We have also analysed the velocity field dependence on di- 
rections defined over the surrounding large scale mater dis- 
tribution. For this purpose we estimate the shape tensor of 
a spherical large scale region around dark matter haloes, ex- 
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Figure 2. Mean infall velocity as a function of distance to the central halo for different shape of halo. Upper panels correspond to high 
mass haloes and lower panels to low masses as is indicated in the figure. Left panels show results for prolate haloes and right panels for 
oblate haloes. Different colours indicate infall directions at different angles from the halo principal axes a and c. Red lines in the left 
panels indicate 8 a < n/6, cyan lines 7r/6< 9 a < tt/3 and black lines n/3< 8 a < n/2. Black lines in right panels indicate 9 C < 7r/6, cyan 
lines 7r/6< 9 C < 7r/3 and red lines n/3< 9 C < tt/2. 



eluding the inner region containing the central halo. In the 
same way as for dark matter haloes, the surrounding large 
scale distribution can be characterised by triaxial ellipsoids 
with positive or negative triaxiality parameters. It can be 
said that LSSR with P > (P < 0) correspond to "fil- 
amentary" ("sheet-like") matter distributions. We use this 
surrounding large-scale structure at the present time (z = 0) 
to characterise the large-scale environment. We expect the 
presence of correlations between the LSSR distribution, and 
the surrounding velocity field at the present time. We stress 
the fact that the dynamical timescale of this volume is large 
enough to contain information on the accretion process of 
the central halo across the halo evolutionary stages. Thus, 
compared to the halo crossing time, the LSSR inertia prin- 
cipal axes can be taken as quasi stationary. For instance, 
a typical halo of 10 13 Mq has a velocity dispersion around 
300 km/s, with a crossing time at virial radius of 1.4 Gyrs, 
whereas a particle with typical velocity of 400 km/s takes 
75 Gyrs to reach a position 10 Mpc h" 1 distant. 

As described above, we characterise the LSSR by means 
of the shape tensor as defined in equation \T\ For the cal- 
culation we use particles between 2 and 15 Mpc h~ from 
halo centres, with a minimum separation threshold so as to 



avoid the central halo particles. In order to fix the maximum 
distance of this spherical shell we have considered several 
thresholds and finally selected the specific value that shows 
the stronger signal. Nevertheless we obtain similar results for 
the different distances considered. The principal axes are la- 
belled with A, B and C, for the major, intermediate and 
minor axes, respectively. We use capital letters in order to 
avoid confusion with the axes a, b and c, which are used to 
characterise the intrinsic halo shapes. Given that the LSSR 
shape tensor represents the second moment of the matter 
distribution around a given halo centre, the A axis corre- 
sponds to the direction along the densest direction, whereas 
the C axis corresponds to the lowest density regions. 

As in the previous section (see equation [2]), we de- 
fine a triaxiality parameter for the surrounding distribu- 
tion using the eigen-values of the LSSR shape tensor, P = 
log((C / B)/(B /A)). In the top panels of Figure Q] we show 
the distribution of P and C/A (dashed black lines). The 
cases with positive parameter, P > 0, (i.e. C/B » B/A, 
and then C < B « A) correspond to structures resem- 
bling to some extent ellipsoidal prolate distributions, here- 
after we shall refer to such regions simply as prolate LSSR; 
filamentary structures are included in this category. On 
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Figure 3. Mean infall velocity as a function of distance to the central halo for different LSSR. The upper panels correspond to high 
mass haloes and the lower panels to low masses as is indicated in the figure. Left panels show results for prolate-LSSR and right panels 
for oblate-LSSR. Different colours indicate infall directions at different angles from the LSSR principal axes A and C . Red lines in the left 
panels indicate 9a < t/6, cyan lines 7r/6< 9a < t/3 and black lines 7r/3< 6a < vr/2. Black lines in the right panels indicate 9q < t/6, 
cyan lines 7r/6< 9c < f/3 an d re d lines 7r/3< 9q < tt/2. 



the other hand, the condition P < (i.e. B/A » C/B, 
so C << B < A ) selects LSSR regions which resemble 
oblate structures. Hereafter we shall call such regions sim- 
ply oblate LSSR. Structures such as pancakes or walls are 
examples of this class. We acknowledge that through this 
analysis of triaxiality parameter we cannot obtain a com- 
plete topological description of the LSSR as would be ob- 
tained, for example, using a Minkowski functional analysis. 
However, given the statistical scope of the current analy- 
sis, we restrict our study to regions with large ellipticities 
(C/A < 0.7), in order to diminish under-determinations of 
the corresponding characteristic axes A and C for both pro- 
late and oblate LSSR samples. In order to assess the align- 
ment between the LSSR distribution and the halo shapes 
we compute the angle 9 between their major axes; the re- 
sult is shown in red lines in the lower panel of Figure [T] 
Similarly the same panel shows the angles between the 
minor axis of halos and LSSR (black lines). As can be 
seen, the major and minor axes of the halo shape and the 
LSSR are high ly aligned, in qualita t ive agreement with pre- 



vious results (fColberg et al.1 2005b jKasun. &i Evrardl 



iBasilakos etaf. 2006; Allgood et al. 2006; Altav et al.1 



2005; 



2006; 



Araeon-Calvo et al.l 120071 : iBrunino et aD 120071 : iBett et al.1 
20071 ; IZhang et al.ll2009l L 



Using a similar analysis to that performed in the pre- 
vious subsection we compute averages of the radial velocity 
field as a function of the distance to the halo centres for 
different directions. The direction is selected using the angle 
(9) subtended by the halo-to-particle position vector with 
respect to the axes A or C, for prolate or oblate LSSR, re- 
spectively. The results are shown in Figure [3] Depending on 
the triaxiality parameter of the LSSR region, the mean ve- 
locity is computed taking into account particles with angles 
9 < 7r/6 with respect to the axis A and C, red (black) lines 
in left (right) panels of Figure [3] On the other hand, we also 
compute the mean radial velocity field in the perpendicular 
directions to these axes requiring 7r/3 < 9 < tt/2, black (red) 
lines in left (right) panels of Figure[3] Intermediate positions 
for both, prolate and oblate LSSR shapes, are obtained in 
the interval 7r/6 < 9 < n/3 (cyan lines). 

The left and right panels in Figure [3] correspond to pro- 
late and oblate LSSR regions, respectively, whereas the up- 
per and lower panels show results for low and high mass 
haloes, respectively. For simplicity we only show the results 
corresponding to the lowest and highest mass samples. As 
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can be seen in the left panels of Figure [3] (prolate envi- 
ronments) particles perpendicular to the major axis (black 
lines) exhibit noticeably higher velocities than particles par- 
allel to this axis (red lines). A similar behavior can be seen 
in the right panels on Figure [3] (oblate LSSR), where parti- 
cles tend to fall faster from directions parallel to the minor 
axis C. Both results are consistent and indicate the prefer- 
ence for a small velocity infall along the characteristic axis 
of both, oblate and prolate surrounding structure. In other 
words, particles show a faster infall when they come from 
low density regions which can be related to the direction of 
the smaller shape tensor eigenvalues. The mean infall ve- 
locity reaches a maximum of approximately 300 km/s and 
exhibits values higher than 400 km/s for the more massive 
haloes. As the distance from halo centre increases, the in- 
fall velocity monotonously decreases. Along the major axis 
direction, particles arrive slower with radial velocities that 
approach constant values or even increase at greater dis- 
tances. 

3.3 Local density and the LSSR 

It is widely accepted that mass follows anisotropic 
streaming mot ions toward cluste rs , preferentially via 
the filaments dColberg et alJ Il999l: (Splinter et all 1 19971 : 
Ivan Haarlem fc van de WeygaertT 1993b). Also, some au- 
thors claim that higher density particles exhibit higher 
infall velocities and coherence than particles in lower 
density environments (se e for instance IPivato et al.l 120061 ; 
IGonzalez fc Padillall2009l ). This behaviour seems to be op- 
posite to the correlation observed in the previous subsection, 
where we found larger infall velocities along the minor axes 
for both, haloes and the surrounding LSSR. Therefore, we 
now explore the dependence on the local density field of the 
mean radial velocities and their anisotropy. 

We define a local density parameter (p) for each parti- 
cle, which measures the mass density in a sphere of radius 
equal to the distance to the 100th nearest neighbour, centred 
around the particle. We also define a mean local density p 
as the average of density parameter logarithms (Log(p)), its 
value corresponds to 5 particles/Mpc 3 for this simulation. 

We separate the particles in the simulation according 
to their associated local density, into high and low density 
samples, defined by p > p and p < p, respectively. The 
left panel of Figure [4] shows the mean radial velocity as a 
function of the distance to halo centres, for particles in low 
and high density regions. We also analyse the dependence 
of these results on the LSSR axis directions. For simplicity 
we only show results for prolate-LSSR. 

The velocity curves shown in this figure correspond 
to low mass central haloes (3.5 x lO 12 M0h _1 < rrihaio < 
4.5 x 10 1 MqIi -1 ), solid lines represent results obtained by 
using particles in dense local environments (p > p), and 
dashed lines correspond to particles in underdense local re- 
gions (p < p). The different colours in Figure [4] correspond 
to different LSSR directions; red lines show the mean veloc- 
ity for particles close to the major LSSR axis (6 a < tt/6) 
whereas black lines show the mean velocity for particles with 
tt/3 < 6 A < tt/2. 

By inspection of the left panel of Figure^ it can be seen 
for both, parallel and perpendicular directions to the LSSR 
major axis (red and black lines respectively), that high den- 



sity particles (solid lines) exhibit higher velocities than low 
density particles (dashed lines). In addition, regardless of 
the particle local density, lower radial velocities are obtained 
along directions parallel to the LSSR major axis. Therefore 
we recover the results obtained in the previous subsection re- 
garding the general dependence of the radial velocity field on 
the LSSR. This figure also reveals a significant dependence 
of radial velocities on the local density confirming previ- 
ous c ommunications (|Pivato et al ]|2006l ; IGonzalez fc Pa dilla 
l2009h . 

Galaxies and groups are expected to flow onto large 
clusters. Thus, it is interesting to study the mean radial 
velocity of haloes in the surroundings of the most massive 
systems. For the purpose of the analysis of halo dynamics 
we have considered similar conditions than those applied to 
particles. 

For haloes we redefine the parameter characterising the 
local environment; the halo local density parameter is calcu- 
lated using the 100th nearest particle neighbour. We exclude 
particles within 3 virial radii in order to avoid the internal 
regions. 

We consider haloes in prolate LSSR regions and com- 
pute the averaged radial velocities of high and low local 
density haloes in preferred LSSR directions. The results are 
shown in the right panel of Figure [3] and consist on the ra- 
dial velocity of low mass haloes (7.2x lO 11 M0h^ 1 < rrihaio < 
4.5 x 10 12 MQh~ 1 ) in the surroundings of high mass haloes 
(1.9 x lO 13 M h _1 < m ha i ). Solid (dashed) lines symbolise 
the infall of haloes in high (low) local densities while red and 
black lines correspond to haloes along, and in the direction 
perpendicular to, the major axis. From a general point of 
view we obtain radial velocities with similar characteristics 
than the velocity field traced by mass particles, namely, low 
density haloes show higher velocities than haloes in dense 
regions. We see an important distinction between infall of 
material in the different LSSR directions, where velocities 
along the major LSSR axis direction are larger than along 
the perpendicular direction. 

From the comparison of both panels of Figure [3] it can 
be seen that the most prominent difference resides on the ve- 
locity amplitudes. Halo velocities are systematically higher 
than for particles. This difference is mainly be due to the 
higher mass of the central haloes, and it is consistent with 
the differences detected in Figure [2] 

We have also explored the infall velocities around cen- 
tral haloes of different masses, and for different LSSR 
shapes, and find consistent results (including higher veloci- 
ties for higher masses). 



4 VELOCITY DISPERSIONS AND 
LAMINARITY OF THE FLOW 

In the previous sections we reported a significant dependence 
of the peculiar velocity field on LSSR directions and local 
density. Our results indicate that radial peculiar velocities 
deliver information on the large scale anisotropies suggesting 
that the velocity dispersions could also be correlated. 

We are now interested in analysing the coherence of 
motions around haloes and its relation with the infall of 
mass onto haloes. To do this analysis we define a laminar- 
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Figure 4. Left panel: Mean infall velocity as a function of distance for particles around low mass haloes (3.5 X IO^MqIi" 1 < mhalo < 
4.5 X 10 12 MQh _1 ) and prolate-LSSR (C/A < 0.7 and P > 0.0) . Solid lines correspond to local densities higher the mean local 
density and dashed lines correspond to local densities lower than mean, as is indicated in the figure. Red lines indicate 9 a < t/6, 
whereas black lines indicate 7r/3< 8\ < 7r/2. Right panel: Mean radial velocity as a function of the distance for low mass haloes 
(7.2 X 10 11 M Q h- 1 < m halo < 4.5 X lO^Meh" 1 ) around high mass haloes (1.9 X 10 13 M Q h- 1 < m halo ) and prolate-LSSR (C/A < 0.7 
and P > 0.0). Lines and colours are as in the left panel. 




Figure 5. Left panel: Mean laminarity parameter (the ratio between the radial velocity and the velocity dispersion in the radial 
direction) as a function of distance for particles around low mass haloes (3.5 X 10 12 Mq1i -1 < m ha i < 4.5 X 10 12 Mq1i — 1 ) and prolate- 
LSSR {C/A < 0.7 and P > 0.0). Solid lines correspond to high local densities and dashed lines to low local densities, as is indicated in the 
figure. Red lines correspond to 8a < t/6, whereas black lines correspond to tt/3< 8a < tt/2. Right panel: Mean laminarity parameter 
as a function of distance for particles around high mass haloes (m^a^ < 1.9 X 10 13 Mq1i -1 ) and prolate-LSSR, lines and colors are as in 
the left panel. 



8 L. Ceccarelli et al. 



I _ _! J _'-L 



prolate-LSSR 




~- 0.6 



0.4 '- 



O<0 A <3O m Ml =4xl0 18 M ffl h-> 3" 
6O<0 A <9O cenler 



_L 



J L 



_L 



J L 



O<0 A <3O m =4xl0 18 M a h- 1 
6O<0 A <9O center 
... I .... I .... I 



5 10 
r/Mpc h _1 



15 



5 10 
r/Mpc h -1 



15 



Figure 6. Upper panels: Velocity dispersions as a function of distance for particles of high and low local densities (left and right panels 
respectively) around low mass halo centres (3.5 X 10 12 Moh~ 1 < mi a j < 4.5 X 10 12 Moh~ x ) and prolate-LSSR. Solid and dashed lines 
indicate radial and perpendicular velocity dispersion respectively, and red lines correspond to 6a < t/6, whereas black lines correspond 
to 7r/3< 8a < 7r/2. Low panels: Ratio between radial and total velocity dispersions as a function of distance for particles around low 
mass haloes and prolate-LSSR. Left (right) panels correspond to high (low) local densities and colours indicate position angles as in the 

upper panels. In the low panels the dotted lines show a \\l + a \ = 1/V2 



ity parameter as the ratio between radial velocity and its 
dispersion. 

In order to obtain a suitable laminarity parameter we 
compute the local velocity dispersion for each halo using 
all the particles (haloes) within spherical shells. As it was 
described above we consider separately the angle subtended 
by the direction to a particle or halo and the principal axes 
of the LSSR. We also take into account the local density. 

In Figures [5] and we show the laminarity parameter 
as a function of the distance to the central halo, for mass 
particles and small haloes respectively. We only show the 
results obtained for central haloes in a prolate LSSR. In 
Figure [S] we select low and high mass central haloes (left 
and right panels respectively), whereas in Figure [7] we show 
velocities of low mass haloes around massive haloes. 

As can be seen in the left panel of Figure [5] along per- 
pendicular directions to the LSSR major axis (black lines) 
particles in either high or low local densities (solid and 
dashed lines) present similar laminarity parameters. In the 
case of particles along the direction parallel to the LSSR ma- 
jor axis (red lines), high density particles (solid lines) show 
higher laminarity parameters than those at low local densi- 
ties (dashed lines). We also obtain similar results using halo 
velocities instead of particles (Figure [7]). When considering 



particle velocities around high mass haloes (right panel of 
Figure [SJ , it can be noticed that along directions perpen- 
dicular to the LSSR major axis (black lines) low density 
particles (dashed lines) show higher laminarity parameters 
than those at high local densities (solid lines). 

From these results we can infer that the flow along the 
different directions becomes more laminar (i.e. v/a > 1) 
closer to the halo centre. This effect is significantly more 
pronounced for infalling clumps than for infalling individ- 
ual particles. This laminarity is also larger along the plane 
perpendicular to the LSSR major axis indicating that, be- 
sides the fact that this direction shows the largest infall, it 
is also the less turbulent; the material is more smoothly ac- 
creted. The implications of this behaviour for halo structure 
could reside in the dynamical memory of the systems of this 
quasi stationary process which acts on time-scales of halo 
formation. 

We also explore possible differences of the peculiar ve- 
locity field along the directions parallel and perpendicular 
to the halo-centric direction. In order to analyse this we cal- 
culate the 1-D velocity dispersion of galaxies in both the ffii 
and, a± directions, respectively, on the frame of the mean 
infall stream. Our results are shown in Figures [6] and [8] 
where we plot the velocity dispersions (top panels) and the 
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Figure 8. Upper panels: Velocity dispersions as a function of distance for low mass haloes (0.7 X lO 12 M0h — 1 < r7i^ a ; < 4.5 X 
10 12 Mq1i — 1 ) of high and low local densities (left and right panels respectively) around high mass halo centres (m^a^ < 1.9 X 10 13 Mq1i -1 ) 
and prolate-LSSR. Solid and dashed lines indicate radial and perpendicular velocity dispersion respectively, and red lines correspond 
to 8a < 7r/6, whereas black lines correspond to ir/3< 6 a < ir/2. Low panels: Ratio between radial and total velocity dispersions as a 
function of distance for low mass haloes around high mass haloes and prolate-LSSR. The left (right) panel corresponds to high (low) 

local densities and colours indicate position angles as in the high panels. In the low panels the dotted lines show ctm /- / <rjf + = 1/V2. 



ratio between radial and total velocity dispersions (hereafter 
the relative dispersion, bottom panels) as a function of dis- 
tance, for both particles and haloes. Different colours indi- 
cate the position relative to the major axis (red lines for 
particles/haloes along the major axis and black lines for the 
perpendicular directions). Left and right panels correspond 
to high and low local densities, respectively. Different lines 
in the top panels indicate radial and perpendicular velocity 
dispersions (solid and dashed lines respectively). For com- 
parison, the expected value for the relative dispersion in an 

isotropic case is f || / y °f + 17 ± = Vv^j an d is shown as 
horizontal dotted lines in the bottom panels. 

As can be seen from the top panels of Figures [6] and [8] 
the dispersion of radial velocities (solid lines) is smaller than 
that of velocities along the perpendicular direction (dashed 
lines). This behaviour can also be seen from bottom panels 
of Figures [6] and [8] where the relative dispersion is lower 
than l/%/2, regardless of the local density, LSSR shape and 
distance to the central halo. We find this behavior out to 
distances as large as 15 Mpc, indicating that the radial infall 
dominates over random motions. 

From the top panels of Figures [5] and [S] it can be seen 
that the velocity dispersion for high local densities (left pan- 
els) is higher than for low local densities (right panels), this 



behaviour is stronger for particle velocities (Figure [B). With 
respect to the LSSR principal axis, the velocity dispersion 
along the mayor axis (red lines) is higher than in its per- 
pendicular direction (black lines) for high local densities at 
distances larger than 5 Mpc h _1 (left top panels of Figures[5] 
and [8}. In general, the velocity dispersion tends to decrease 
towards smaller distances from the central halo. 

The lower panels of Figures [S] and [5] which correspond 
to low local densities, show that the relative dispersion along 
the major axis decreases as the distance to the central 
halo decreases. At smaller distances from the central halo 
(r < 10 Mpc h _1 ) the relative dispersion along the major 
axis (red lines) is smaller than along the perpendicular di- 
rection (black lines), whereas at larger distances (r < 13 
Mpc h _1 ) we find the opposite behaviour. We note that the 
dependence of the relative dispersion on the direction with 
respect to the LSSR major axis tends to be smaller as the 
distance increases (as expected due to the tendency towards 
homogeneity of the LSSR) . 

The relative dispersion for high local densities (left bot- 
tom panels of Figures [5] and |SJ) shows similar behaviours for 
the different LSSR directions and, as it can be seen, there is 
little dependence on the distance to the central halo with a 

roughly constant value tr^/./a?, + a\ — 0.6. We do not de- 
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Figure 7. Mean laminarity parameter as a function of dis- 
tance for low mass haloes (7.2 X 10 11 Mq1i — 1 < raj a | < 4.5 X 
lO 12 M h _1 ) around high mass haloes (1.9 X lO 13 M h _1 < 
m ha i ) and prolate-LSSR (C/A < 0.7 and P > 0.0). Red lines 
indicate 8a < t/6, whereas black lines indicate ir/3< 9a < vr/2. 
The two pairs of solid and dashed lines correspond to high and 
low local densities as is indicated in the figure. 

tect a strong dependence of the relative dispersion on local 
density and LSSR. A more detailed analysis shows a slight 
tendency of the relative dispersion to be smaller along the 
radial than the perpendicular direction, for both, particles 
and haloes, as is seen in the right bottom panels of Figures 
E]and El 



5 INFALL AND OUTFLOW 

In this section we provide further statistical properties of 
matter infall to, and outflow from haloes. For this aim we 
have considered separately particles with negative (infalling) 
and positive (outflowing) radial velocity components. There- 
fore, these two groups of particles consist in local diverging 
and converging flows onto the halos. 

Velocity moments for each group of converging or di- 
verging particles are obtained in a similar way as the one 
used to obtain the shape tensor (Equation [T]), using the 
components of the peculiar velocity vector of each parti- 
cle relative to the centre of mass velocity of the surrounding 
LSSR (within a sphere of 15 Mpc h _1 from the halo centres) 
instead of the position vector. 

Then, for each group of particles we obtain a triaxial 
ellipsoid that approximates the velocity distribution. The 
corresponding major, intermediate and minor principal ve- 
locity axes are referred to as a, /3 and 7. We analyse the 
orientation of these two velocity ellipsoids (for infalling and 
outfalling particles), with respect to the principal axes of the 
surrounding structure. As in previous sections, we refer to 
these spatial axes as A, B and C, for the major, intermediate 



and minor axes, respectively. We compute the angles sub- 
tended by the principal axes of shape and velocity ellipsoids 
and summarise our results in Figure El This figure shows 
the distributions of the cosine of the angles between LSSR 
and velocity axes, for converging (black lines) and diverging 
particles (red lines). Each panel in this figure corresponds to 
a different pair of axes, the three upper panels correspond 
to angles between the LSSR major axis and the velocity 
major, intermediate and minor axes (panels a, b and c, re- 
spectively). Similarly, the three middle panels correspond to 
the LSSR intermediate axis and the three lower panels for 
the LSSR minor axis, as indicated in the figure. 

As can be seen in the left upper panel of this figure, 
the major LSSR and velocity axis tend to be perpendicu- 
lar for converging particles (black lines in panel (a)), and 
parallel for diverging particles (red lines in panel (a)). The 
right upper panel of this figure shows that the major LSSR 
axis tends to be aligned with the minor velocity axis for 
converging particles (black lines in panel (c)) whereas for 
diverging particles they tend to be perpendicular (red lines 
in panel(c)). Regarding the angles between the LSSR major 
axis and the velocity intermediate axis there are no signifi- 
cant differences between converging and diverging particles 
(black and red lines in panel (b)). By inspection of the lower 
panels of this figure, which correspond to the angles between 
the LSSR minor axis and the velocity ellipsoid, it can be 
seen, consistent with the results in the upper panels, that 
the minor LSSR axis and the major velocity axis are aligned 
for converging particles and perpendicular for diverging par- 
ticles (black and red lines in panel (i)). 

The main results shown in Figure El indicate that LSSR 
major axes tend to be parallel to major (minor) velocity 
axes for diverging (converging) particles. The alignment sig- 
nals for the converging flow is consistent with the results 
derived in the previous sections, namely a preference for 
larger amplitudes of the infalling velocity field along direc- 
tions perpendicular to the principal axes of the surrounding 
structure. 

Regarding the outflowing velocity particles we notice 
the opposite behaviour, the largest amplitude occurs along 
the LSSR major axis direction. Thus, in part, the reduced 
mean infall along the LSSR major axis could be caused by 
the preference of outflowing particles along this direction. 
We can interpret this trend as generated by particles ex- 
pelled from the virialised structures, which could occur more 
easily along the major axis where the gravitational potential 
shows a more gentle decrease than along the minor axis di- 
rection. Thus, a possible implication of our results would be 
the presence of galaxies that passed close to the halo cen- 
tres and are outflowing along the surrounding filamentary 
and planar structures. 



6 CONCLUSIONS 

Using a numerical simulation we have performed a detailed 
analysis of the peculiar velocity field around dark matter 
haloes. We have explored the relation between the flows of 
mass towards haloes and the shapes of the central haloes, 
their large-scale surrounding region (LSSR) shape, and their 
orientation. 

We found a significant anisotropy in the velocity field 
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Figure 9. Histograms showing the vectorial products between principal axes derived from the infall velocity tensor and the LSSR for 
particles (dashed black lines). Each panel represents the angle between different pairs of vectors as is indicated in the corresponding 
panel in the figure. A, B, and C represent the major, intermediate and minor axes of the LSSR, respectively, a, j3 and 7 represent the 
major, intermediate and minor axes of the infall velocity ellipsoid, respectively. The corresponding histograms derived from the outflow 
velocity tensor and the LSSR are shown as solid red lines. 



which correlates with the surrounding LSSR principal axes. 
The amplitude of the velocity field of infalling particles onto 
dark matter haloes is maximum along the direction of under- 
dense surrounding regions (LSSR minor axis), whereas along 
the direction of dense LSSR regions, it exhibits the smallest 
velocities. With respect to the shape ellipsoid of the central 
haloes, we found that the amplitude of the infall velocity 
field along the halo minor axis is larger than that along the 
major axis. We find consistent results for general triaxial 
haloes, and for both prolate and oblate systems. Therefore, 
there is a clear preference for matter to infall towards haloes 
along regions of low "global" density. 

When looking at the dependence with local density, ei- 
ther measured in cubic cells in the simulation volume for in- 
dividual particles, or using the 100th closest particle neigh- 
bour outside 3 virial radii for tracer haloes, we find that 
objects in local high density regions infall faster. 

By defining a suitable laminarity parameter we find that 



the velocity field is less turbulent along the direction to- 
wards underdense surrounding regions, with a ratio between 
mean flow velocity and velocity dispersion of order unity and 
nearly constant up to scales of 15 Mpc h _1 . 

The flow of particles moving away from haloes are found 
to be highly anisotropic. The major axis of the outflowing 
velocity field is aligned with the surrounding LSSR major 
axis. 
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